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A  portable  proton  exchange  membrane  (PEM)  fuel  cell-battery  power  system  that  uses  hydrogen  as  fuel 
has  a  higher  power  density  than  conventional  batteries,  and  it  is  one  of  the  most  promising  environ¬ 
mentally  friendly  small-scale  alternative  energy  sources.  A  general  methodology  of  modeling,  control 
and  building  of  a  proton  exchange  membrane  fuel  cell-battery  system  is  introduced  in  this  study.  A  set 
of  fuel  cell-battery  power  system  models  have  been  developed  and  implemented  in  the  Simulink  envi¬ 
ronment.  This  model  is  able  to  address  the  dynamic  behaviors  of  a  PEM  fuel  cell  stack,  a  boost  DC/DC 
converter  and  a  lithium-ion  battery.  To  control  the  power  system  and  thus  achieve  proper  performance, 
a  set  of  system  controllers,  including  a  PEM  fuel  cell  reactant  supply  controller  and  a  power  manage¬ 
ment  controller,  were  developed  based  on  the  system  model.  A  physical  100W  PEM  fuel  cell-battery 
power  system  with  an  embedded  micro  controller  was  built  to  validate  the  simulation  results  and  to 
demonstrate  this  new  environmentally  friendly  power  source.  Experimental  results  demonstrated  that 
the  100  W  PEM  fuel  cell-battery  power  system  operated  automatically  with  the  varying  load  conditions 
as  a  stable  power  supply.  The  experimental  results  followed  the  basic  trend  of  the  simulation  results. 
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1.  Introduction 

Fuel  cells,  a  promising  next-generation  power  source,  are  widely 
used  in  both  automotive  and  stationary  applications  due  to  their 
high  power  density  and  low  emissions.  A  proton  exchange  mem¬ 
brane  (PEM)  fuel  cell  is  an  electrochemical  device  in  which  the 
energy  of  an  electrochemical  reaction  is  converted  directly  into 
electricity  by  combining  hydrogen  fuel  with  oxygen  from  the  air 
with  heat  as  the  bypass  [1  ]. 

A  PEM  fuel  cell  that  operates  in  low  temperature  conditions  is 
suitable  for  the  portable  power  application.  However,  PEM  fuel  cells 
have  issues  working  independently  as  a  power  system.  A  rapid  volt¬ 
age  drop  and  a  slow  response  to  the  load  demand  are  the  primary 
critical  issues,  as  most  electronic  devices  require  stabilized  power 
and  fast  transient  response.  Thus,  energy  storage  devices,  such  as 
fast  charge  and  discharge  battery  and  power  regulation  devices, 
such  as  DC/DC  converters,  are  required  to  work  with  the  PEM  fuel 
cell  to  form  a  power  system  to  provide  regulated,  fast-response 
power  to  a  variable  load. 

A  typical  PEM  fuel  cell-battery  power  system  consists  of  a  PEM 
fuel  cell  stack  as  the  main  power  source,  a  buck  or  boost  DC/DC 
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converter  as  the  power  conditioning  device  and  a  secondary  battery 
power  source. 

For  the  best  understanding  of  the  PEM  fuel  cell  operation,  PEM 
fuel  cell  modeling  is  required  to  numerically  examine  the  PEM 
fuel  cell  system.  The  basic  concept  of  the  fuel  cell  and  the  general 
mathematical  model  has  been  introduced  in  the  literature  [1-3]. 
Previously,  many  improved  PEM  fuel  cell  mathematical  models 
have  been  published  [4-12].  Most  of  these  are  intended  for  PEM 
fuel  cell  steady-state  simulation  or  for  the  sizing  of  PEM  fuel  cell 
parameters.  Due  to  the  complexity  and  heavy  calculative  load  of 
these  models,  they  are  not  suitable  for  control  design  pursuers.  To 
control  and  operate  a  PEM  fuel  cell  system  at  an  optimal  condition 
and  to  maintain  high  performance  of  the  fuel  cell  system,  a  sim¬ 
plified,  control-oriented  model  is  required  for  the  control  design. 
Pukrushpan  et  al.  first  developed  a  dynamic  model  that  is  suit¬ 
able  for  a  control  study  of  fuel  cell  systems  [13].  They  simplified 
existing  models  to  make  a  sufficiently  simple  model  that  can  simul¬ 
taneously  address  the  major  PEM  fuel  cell  transient  and  behavior 
and  significantly  reduce  the  complexity  and  computational  load. 
This  simplified  model  ensured  that  the  controls  for  the  PEM  fuel 
cell  could  be  rapidly  developed.  Later,  several  improved  control- 
oriented  models  were  published  [14-16].  These  models  use  the 
same  principle  to  derive  the  PEM  fuel  cell  model  and  add  either 
a  temperature  effect  model  or  a  humidifier  model  to  the  control- 
oriented  PEM  fuel  cell  model. 

After  the  control-oriented  PEM  fuel  cell  model  was  well  devel¬ 
oped  and  widely  implemented,  the  control  design  for  PEM  fuel  cell 
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Nomenclature 

a,  water  activity 

Afc  fuel  cell  active  area  (cm2) 

cv  water  concentration  ( mol  m-3 ) 

Cp  specific  heat  capacities  (J  kg-1  K-1 ) 

Do  diffusion  coefficient  (m2  s-1 ) 

F0  open  circuit  voltage  (V) 

F  Faraday  constant  (F=  96,485) 

i  current  density  (A  cm-2 ) 

/0  exchange  current  density  (A  cm-2 ) 

ibatt  battery  current  (A) 

id  low  frequency  current  dynamics  (A) 

imax  limiting  current  (A) 

/  current  (A) 

m  mass  (kg) 

m  mass  flow  rate  (kg  s-1 ) 

M  molecular  mass  (kg  mol-1 ) 

n  number  of  cells 

rid  electro-osmotic  coefficient 

p  pressure  (Pa) 

Psat  vapor  saturation  pressure  (Pa) 

Q.  maximum  battery  capacity  (Ah) 

Qext  extracted  capacity  (Ah) 

Qexp  exponential  capacity  ((Ah)-1 ) 

R  gas  constant  (=8.314J  mol-1  K-1 ) 

Rc  equivalent  series  resistor  of  the  capacitor  (£2) 

Ri  body  resistor  of  the  inductor  ( £2 ) 

tm  thickness  of  membrane  (£2) 

T  temperature  (K) 

V  over  potential  (V) 

Vexp  exponential  voltage  (V) 

W  work  (J) 

crm  membrane  conductivity  (£2-1  cm-1 ) 
y  ratio  of  the  specific  heat  capacities  of  the  gas 

pCp  isentropic  efficiency 

Xo2  oxygen  excess  ratio 

Ah2  hydrogen  excess  ratio 

Am  membrane  water  content 

0  relative  humidity 

&>atm  air  humidity  ratio 

Subscripts 
a  air 

act  activation  loss 

atm  atmosphere 

an  anode 

ca  cathode 

cell  single  cell 

cone  concentration  loss 

fc  fuel  cell 

gen  generated 

H2  hydrogen 

in  inlet 

membrane  across  membrane 
N2  nitrogen 

ohm  ohmic  loss 

out  outlet 

02  oxygen 

reacted  reacted 
st  stack 

v  water  vapor 


system  was  simplified.  All  of  the  known  control  methods  can  be 
applied  to  PEM  fuel  cell  applications.  A  number  of  publications 
have  been  released  [13,15,17,18]  with  different  control  strategies 
(e.g.,  linear,  non-linear,  and  fuzzy-logic)  or  for  different  control  pur¬ 
poses  in  fuel  cell  operation  (e.g.,  flow  control,  humidity  control,  and 
pressure  control). 

As  the  major  application  of  the  PEM  fuel  cell  stack,  the  PEM  fuel 
cell-battery  power  system  is  widely  used  in  automotive  or  station¬ 
ary  applications.  The  PEM  fuel  cell-battery  power  system,  including 
the  system  design,  the  modeling  and  the  power  management  sys¬ 
tem,  has  been  discussed  in  the  literature  [19-24].  Elowever,  the 
fuel  cell  dynamics  flow  behaviors  were  typically  neglected  and  a 
PEM  fuel  cell  equivalent  circuit  was  used  to  represent  the  PEM 
fuel  cell  model  when  the  PEM  fuel  cell  stack  was  modeled.  This 
could  reduce  the  complexity  of  the  model;  however,  the  dynamic 
flow  behavior  of  the  fuel  cell  is  the  key  component  in  PEM  fuel 
cell  performance.  Neglecting  the  dynamic  flow  behavior  of  the  fuel 
cell  could  result  in  a  significant  reduction  in  the  model  accuracy. 
The  PEM  fuel  cell-battery  power  systems  in  these  publications 
may  have  different  topologies;  however,  the  power  management 
strategies  all  have  a  similar  philosophy,  which  is  a  balanced  power 
distribution  between  the  fuel  cell  and  the  battery  to  satisfy  the  load 
demand. 


2.  PEM  fuel  cell  stack  model 

Many  PEM  fuel  cell  stack  models  have  been  previously  proposed. 
Currently,  the  best  control-oriented  PEM  fuel  cell  stack  model 
was  developed  by  Pukrushpan  et  al.  They  developed  a  simplified 
dynamic  model  that  is  suitable  for  a  fuel  cell  control  study  [13]. 
In  the  model,  the  temperature  of  the  fuel  cell  stack  was  assumed 
to  be  constant.  This  assumption  reduced  the  order  of  the  PEM  fuel 
cell  stack  model,  the  simulation  time  and  the  complexity.  In  con¬ 
trast,  the  temperature  of  the  fuel  cell  stack  can  be  controlled  using 
a  separate  thermal  management  system.  Thus,  this  assumption  is 
reasonable. 

The  PEM  fuel  cell  stack  model  introduced  in  this  study  is  pri¬ 
marily  based  on  this  model.  A  PEM  fuel  cell  stack  model  consists 
of  two  components:  the  PEM  fuel  cell  thermodynamics  and  the 
electrochemical  reaction. 


2.1.  Dynamic  flow  model  in  PEM  fuel  cell 

For  the  operation  of  the  PEM  fuel  cell,  the  thermodynamic  model 
is  the  major  factor  that  dominates  the  fuel  cell  performance  and 
stability. 

In  this  model,  the  law  of  mass  conservation  was  used  to  address 
the  mass  flow  rate  of  oxygen,  nitrogen,  and  vapor  in  the  cathode 
and  the  hydrogen  and  vapor  in  the  anode  [13]. 


2.1.1.  Mass  flow  governing  equations 

Mass  conservation  law  was  used  to  derive  the  mass  flow  rate  of 
oxygen,  nitrogen,  and  vapor  in  the  cathode  and  the  hydrogen  and 
vapor  in  the  anode  [13]. 

On  the  cathode  side,  this  results  in  the  following  equations: 


dmo2Ca  . 

^  —  ^02,ca,in  —  ^02,ca,out  —  m02,  reacted 

^mN2,ca  ^  ^ 

^  —  mN2,ca,in  —  mN2,ca,out 

dmv,ca  ^  ^  [  [-i 

^  —  ^v,ca,in  —  ^v,ca,gen  +  ^v,memberane 


(1) 

(2) 

(3) 


On  the  anode  side,  this  results  in  the  following  equations: 


dm 


H2,an 


dt 


=  ™h2 


l  ^H2,ca,out  +  ^H2,  reacted 


(4) 
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dmy 


dt 


—  ^v,an,in  ^v,  an,  out  +  ^v,  membrane 


(5) 


The  electrochemical  principle  was  used  to  calculate  the  rates  of 
oxygen  consumption,  water  production,  and  the  hydrogen  con¬ 
sumption  rate  from  the  stack  current,  Jst: 


nl  st 

m°2,  reacted  =  ^°2  X  4 p 

(6) 

n/r 

lPlv,c a, gen  =  x  -^r 

(7) 

nl  st 

mH2, reacted  —  ^H2  x 

(8) 

where  n  is  the  number  of  cells  in  the  stack  and  F  is  the  Faraday 
number,  which  is  equal  to  96,485  C  mol-1 . 

2.1.2.  Stack  pressure  change 

In  this  model,  all  gases  were  assumed  to  obey  the  ideal  gas  law 
and  the  temperatures  inside  the  cathode  and  anode  were  assumed 
to  be  equal  to  the  stack  temperature  [13].  Then  the  partial  pressures 
of  oxygen  (po2 ).  nitrogen  (pN2 ),  the  cathode  vapor  (pv,ca),  hydrogen 
(Ph2),  the  anode  vapor  (pv,an),  and  the  relative  humidity  (0)  were 
calculated  using  the  ideal  gas  law  based  on  the  masses  of  oxygen 
(m02),  nitrogen  (mN2),  the  cathode  vapor  (mv,ca),  hydrogen  (mH2), 
the  anode  vapor  (mv,an).  and  the  stack  temperature  (Tfc). 

In  the  cathode,  the  following  equations  were  used. 


Oxygen  partial  pressure: 

Po2,ca  — 


nio2,ca^02^st 


Vca 

Nitrogen  partial  pressure: 

n  niN2)CaFN2Tst 

PN2,ca  -  n 

vca 

Vapor  partial  pressure: 


(9) 


Pv,ca  = 


nTv,caPv^st 

VZ 


Partial  pressure  of  dry  air  in  cathode: 

Pa,ca  =  Po2,ca  +  PN2,ca 

Total  cathode  pressure: 

Pea  =  Pa,ca  +  Pv,c  a 

In  the  anode,  the  following  equations  were  used. 


(11) 

(12) 

(13) 


Hydrogen  partial  pressure: 

^H2,an^H2^st 

PH2'an  “  ^an 

(14) 

Vapor  partial  pressure: 

n7v,anFvTst 

Pv,an  —  w 

van 

(15) 

Total  anode  pressure: 

Pan  —  PH2,an  +Pv,an 

(16) 

2 A. 3.  Relative  humidity  (RH)  inside  stack 
Relative  humidity  in  cathode: 


0ca  = 


Pv,c  a 


Psat(Tst) 

Relative  humidity  in  anode: 

Pv,an 


(17) 


0an  — 


where  psa t  is  the  vapor  saturation  pressure  that  varied  with  tem¬ 
perature.  In  this  model,  Tst  was  a  constant  value;  therefore  the 
vapor  saturation  pressure  was  also  a  constant. 

2 A  A.  Water  transfer  across  the  membrane 

The  water  transfer  behavior  across  the  membrane  is  a  result  of 
two  major  phenomena:  the  electro-osmotic  drag  from  the  anode 
to  the  cathode  and  back-diffusion  from  the  cathode  to  the  anode. 
Here,  the  water  concentration  (cv)  was  assumed  to  change  linearly. 
The  following  equation  was  used  to  capture  this  phenomenon: 

(Cy,av  —  Cv,an)  \ 


^v,  membrane  —  A/fv71fcn 


n  ~  D(> 

F 


(19) 


where  tm  is  the  thickness  of  the  membrane,  nd  is  the  electro- 
osmotic  coefficient,  D M  is  the  diffusion  coefficient,  which  varied 
with  the  water  content  in  the  membrane,  Mv  is  the  vapor  molar 
mass,/\fc  is  the  fuel  cell  active  area,  and  n  is  the  number  of  fuel  cells 
in  the  stack  [13]. 

The  membrane  water  content  (Am)  was  calculated  using  the 
following  equation: 

,  f  0.043  +  17.81a, -39.85a2  +  36.0a3  (0<a,<l) 

*m=  \  14  +  1.4(0, -1)  (0  <  a,  <  3)  (20) 

where  a,  is  the  water  activity,  which  can  be  calculated  using  Eq. 

(21). 


a#  = 


aan  +  tlc  a 


(21) 


The  electro-osmotic  and  diffusion  coefficients  were  calculated 
using  Eqs.  (22)  and  (23): 

nd  =  0. 0029A.ii,  +  0.05Am  -  3.4  x  1(T19 


(10)  =  exp  (24,6  (3I3-J-)) 


where 


Dx  =  < 


10-6  (Am  <  2) 

10“6(1  +2(Am-2))  2  <  Am  <  3 
10_6(— 1 .67(Am  —  3)  3  <  Am  <  4.5 


(22) 

(23) 

(24) 


1.25  x  10“ 


4.5  <  An 


Tfc  is  the  temperature  of  the  fuel  cell. 


2.2.  Electrochemical  model 


The  electrochemical  model  was  primarily  used  to  address  the 
electric  property  of  the  PEM  fuel  cell  stack. 

The  fuel  cell  voltage  model  was  developed  by  subtracting  the 
three  major  over  potentials  from  the  open  circuit  voltage: 


VCell  —  Eq  —  Vact  —  ^ohm  — 


(25) 


where  F0  is  the  cell  open  circuit  voltage,  also  known  as  the  thermo¬ 
dynamic  potential,  was  calculated  using  Eq.  (26). 

E0  =  1.229-0.85  x  10-3(Tfc  -  298.15) 


+  4.3085  x  10“5rfc 


ln(PH2)+  ^ln(Po2) 


(26) 


The  activation  over  potential  was  caused  by  the  slow  charge  trans¬ 
fer  reaction  at  the  surface  of  the  electrodes.  Some  of  the  electrode 
potential  is  used  to  drive  the  electron  transfer  to  match  the  current 
demand;  thus,  the  potential  at  the  fuel  cell  terminal  is  reduced. 

The  activation  over  potential  was  determined  using  the  follow¬ 
ing  equation: 


Psat(Tst) 


(18)  Vact  =  a  In 


(i) 


(27) 
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where  i  is  the  current  density,  i0  is  the  exchange  current  density.  i0 
and  a  are  all  constants  and  they  were  determined  empirically. 

The  ohmic  over  potential  occurs  due  to  the  resistance  to  the 
flow  of  ions  in  the  electrolyte;  the  resistance  to  the  flow  of  elec¬ 
trons  through  the  electrodes  and  the  contact  resistance  at  the  cell 
terminals.  The  ohmic  over  potential  ( Vohm )  can  be  determined  using 
the  following  equation: 


The  inlet  and  outlet  temperatures  and  pressures  are  defined  as 
Ti ,  T2,  pi  and  p2 ,  respectively. 

Mgr* 

where  y  is  the  ratio  of  the  specific  heat  capacities  of  the  gas  (cp/cv). 
The  following  assumptions  were  used  to  simplify  the  calculation. 


Vohm  =  «ohm  =  (28) 

Om 

where  Rohm  is  the  internal  electrical  resistance,  tm  is  the  thickness 
of  the  membrane,  and  the  membrane  conductivity  is  crm,  which  can 
be  calculated  using  the  membrane  water  content  Am  and  the  fuel 
cell  temperature  Tfc. 

The  value  of  the  membrane  water  content  (Am)  varies  between 
0  and  14,  which  is  equivalent  to  the  relative  humidity  from  0%  to 
100%,  respectively. 

<Jm  =  (bnA.m  -bu)  exp  (b2  -  3-))  (29) 

The  constants  bn,  b\2  and  b2  are  typically  determined  empirically 
[13]. 

The  concentration  over  potential  occurs  due  to  a  decrease  in 
the  concentration  of  the  reactants  at  the  electrode-electrolyte 
interface.  A  steady  supply  of  the  reactants  is  required  at  the 
electrode-electrolyte  interface  to  maintain  the  flow  of  electric  cur¬ 
rent.  Due  to  diffusion  or  convection  issues  in  the  electrolyte,  the 
concentration  of  the  reactants  is  not  maintained  at  the  initial  level. 
This  over  potential  ( VConc)  was  calculated  using  the  following  equa¬ 
tion: 

Vconc  =  l  fc274-)c3  (30) 

V  ‘max  / 

where  c2,  C3,  and  imax  are  constants  that  were  determined  empiri¬ 
cally.  The  imax  is  the  limiting  current. 

Because  the  individual  cells  were  stacked  in  series,  the  stack 
voltage  was  calculated  using  Eq.  (31 ): 

Vst  =  n  x  VfC  (31) 

where  n  is  the  number  of  cells. 


2.3.  Air  and  hydrogen  supply  model 


In  this  study,  the  air  was  assumed  to  be  instantaneously  supplied 
to  the  fuel  cell  stack  cathode  by  an  electronic  blower,  and  the  flow 
rate  was  calculated  as  a  function  of  the  required  fuel  cell  current: 


^ca,in  —  (1  +  ^atm)fhdryair  —  (1  +  ^atm) 

*o2 


Xq2Mq 


nht 

2  4F 


(32) 


where  Ao2  is  the  oxygen  excess  ratio  and  co atm  is  the  air  humidity 
ratio. 

The  hydrogen  was  supplied  from  a  hydrogen  tank  and  controlled 
using  a  solenoid  valve.  The  flow  rate  was  also  calculated  as  a  func¬ 
tion  of  the  required  fuel  cell  current: 


m 


an,  in 


=  ^h2Mh 


st 

2  2F 


(33) 


where  Ah2  is  the  hydrogen  excess  ratio. 


3.  Fuel  cell-battery  power  system  auxiliary  model 


3.1.  Air  blower  model 


The  blower  model  introduced  in  this  study  is  a  general  com¬ 
pressor/blower  model  previously  discussed  in  [2]  that  uses  the 
pressure  difference  between  the  blower  inlet  and  outlet  and  the 
blower  efficiency  to  calculate  the  power  consumed  by  the  blower. 


1.  The  heat  generated  by  the  mechanism  was  neglected. 

2.  The  change  in  kinetic  energy  between  the  inlet  gas  and  the  outlet 
gas  were  negligible. 

3.  The  gas  specific  heat  at  constant  pressure  (cp)  was  considered  to 
be  constant  during  the  compression  process. 


Under  the  above  assumptions,  enthalpy  is  the  only  state  that  is 
changed  by  the  mechanical  work,  resulting  in  the  following  equa¬ 
tion: 


W  —  Cp(I2  —  T\  )/TZgas 


(35) 


rhgas  is  the  mass  rate  of  the  compressed  gas.  The  above  formula 
represents  the  isentropic  process  because  the  exit  temperature  for 
the  real  work  will  be  higher  than  the  isentropic  temperature: 

W*=CP(T*-T1)mgas  (36) 

W*  and  T*  are  the  real  work  of  compression  and  the  real  exit  tem¬ 
perature,  respectively.  The  ratio  between  the  isentropic  work  and 
the  real  wok  can  be  calculated;  this  ratio  is  known  as  the  isentropic 
efficiency. 

_  isentropic  work  _  cp(T2  -  h  )mgas  _  T2  -T\ 

^cp  -  real  work  -  cp(T|  -  Ti)mgas  -  T*  -  T\ 

By  combining  Eqs.  (35)  and  (37),  the  equation  shown  below  is 
found. 

h 

Vcp  r*  _  Ti 

The  temperature  difference  is  shown  in  Eq.  (39): 

at  =  Ti  -  r,  =  — 

Vcp 


(pO 


r-Vr 


- 1 


(39) 


y-Vy 


-1 


(38) 


The  combination  of  Eqs.  (35)  and  (39)  results  in  the  following  equa¬ 
tion. 


Power  =  W*  =  cp  — 

hep 


CT’/r 


-1 


hlgas 


(40) 


This  calculated  power  is  the  mechanical  power  required  to  raise 
the  gas  pressure.  To  calculate  the  electric  power  consumed  by  the 
blower,  the  efficiency  between  the  driving  motor  and  the  turbine 
must  be  considered.  In  this  study,  the  efficiency  was  specified  as 
85%,  according  to  the  AMETEK  3.0  BLDC 1 2  VDC  Low- voltage  Blower 
data  sheet. 


3.2.  DC/DC  converter  model 

DC/DC  converters  are  used  to  regulate  the  fuel  cell  output  power 
by  converting  the  rapidly  changing  fuel  cell  output  voltage  to  a  sta¬ 
bilized  value  to  meet  the  requirements  of  electronics.  In  this  study, 
the  selected  fuel  cell  stack  output  voltage  was  less  than  1 0  V.  For  the 
power  to  be  used  for  a  standard  electronic  device,  the  voltage  was 
boosted  to  12  V,  which  is  the  nominal  voltage  for  auto  electronics 
and  small  portable  electronic  devices.  Thus,  the  12  V  boost  DC/DC 
converter  was  chosen  for  this  study,  and  its  efficiency  was  assumed 
to  be  85%  [25]. 

Fig.  1  shows  the  principle  of  a  general  switching  boost  DC/DC 
converter. 
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Fig.  1.  Boost  DC/DC  converter. 


The  ON-OFF  state  demonstrated  for  a  boost  DC/DC  converter  is 
shown  in  Fig.  2. 

The  basic  principle  of  a  boost  converter  includes  2  distinct  states 
(see  Fig.  2). 

In  the  ON-state,  the  switch  (see  Fig.  2)  is  closed,  resulting  in  an 
increase  in  the  inductor  current. 

In  the  OFF-state,  the  switch  is  open  and  the  only  path  offered 
to  the  inductor  current  is  through  the  fly  back  diode,  the  capaci¬ 
tor  (C)  and  the  load  (1?).  This  results  in  the  transfer  of  the  energy 
accumulated  during  the  ON-state  into  the  capacitor. 

Based  on  the  basic  principle  of  this  DC/DC  converter,  a  mathe¬ 
matical  model  was  derived  as  follows  [26]: 


„ dvc  ..  V0 

CW=(  l-uh-j-to 

(41) 

L~df  =v°~^  ~  u^v°  ~  ^L*L 

(42) 

Rvc  RRq  , .  .  . 

V°  ~  R  +  Rc  +  R  +  Rc  ‘L  ‘° 

(43) 

where  RL  is  a  body  resistor  of  the  inductor  and  Rc  is  the  equivalent 
series  resistor  of  the  capacitor. 


3.3.  Battery  model 

The  generic  battery  model  developed  by  Tremblay  [27]  is  a  sim¬ 
ple  battery  model;  however,  it  can  be  used  to  capture  the  primary 
behaviors  of  the  battery.  Thus,  the  lithium-ion  battery  model  intro¬ 
duced  by  Tremblay  [27]  was  employed  in  this  paper. 


On-State 


Off-State 


Discharge  model  (id  >  0): 


/l  (Qext,  id,  ibatt)  -E0-R-  Q  _<^2ext  •  «d  -  R  ■ 

+  l/exp(  —  Qexp  •  Qext)  (44) 

Charge  model  (id  <  0): 


/2(0ext,  id 5  ibatt)  —  Eq  —  R  • 


Q 

Qext  +  0.1  •  Q 


•id— R- 


Q 

Q  —  Qext 


•  Qext 


+  ^exp  •  6Xp(— Qexp  •  Qext) 


(45) 


State  of  charge  (SOC): 


(46) 


where  E0  =  open  circuit  voltage  (V),  R  =  internal  resistance  (£2), 
id  =  low  frequency  current  dynamics  (A),  ibatt  =  battery  current 
(A),  Qext  =  extracted  capacity  (Ah),  Q=  maximum  battery  capacity 
(Ah),  Vexp  =  exponential  voltage  (V),  and  Qexp  =  exponential  capac¬ 
ity  ((Ah)-1). 


4.  System  controls  and  simulation  results 

4A.  Air  supply  control 

Previous  studies  have  proposed  different  types  of  control  strat¬ 
egy  for  the  air  supply  control  in  PEM  fuel  cells.  Most  of  these 
methods  used  different  types  of  feedback  control,  which  collects 
different  system  states  using  many  sensors  and  designs  a  high 
robustness  controller  to  tolerate  system  uncertainties.  In  doing  so, 
the  controller  could  theoretically  produce  accurate  results;  how¬ 
ever,  it  significantly  increases  the  system  cost  and  the  heavy  load 
calculation  is  not  suitable  for  implementation  into  the  embedded 
controllers.  In  contrast,  it  is  known  from  experiments  that  the  most 
serious  current  challenge  for  the  unstable  operation  of  PEM  fuel 
cells  is  a  liquid  water  problem,  which  cannot  be  solved  by  model- 
base  control  alone.  This  is  because  there  is  no  non-CFD  model  that  is 
able  to  address  the  liquid  water  behavior  in  a  PEM  fuel  cell,  and  CFD 
models  cannot  be  used  for  control  development.  Therefore,  in  our 
application,  feed-forward  control  for  the  air  supply  system  was  the 
best  choice.  The  feed  forward  control  in  our  application  presents 
several  advantages. 

i.  It  does  not  introduce  instabilities. 

ii.  It  does  not  rely  on  measurements  of  the  system’s  state,  which 
were  much  slower  than  the  dynamics  of  the  system  in  this  par¬ 
ticular  case. 

iii.  It  relies  on  a  measurement  of  the  system’s  only  disturbance  (Js  t), 
which  we  can  be  obtained  inexpensively  with  high  precision  and 
a  large  bandwidth. 

iv.  It  reduced  the  computational  load  for  the  embedded  controller. 

The  feed  forward  control  algorithm  for  air  supply  control  is 
shown  in  Fig.  3.  The  load  required  current  was  measured  using  a 
current  sensor  as  the  load  following  signal  and  then  converted  into 
a  digital  signal.  The  load  following  signal  acts  as  the  only  input  for 
the  feed  forward  controller.  The  feed  forward  coefficient  was  calcu¬ 
lated  as  a  function  of  the  load  using  Eq.  (32).  After  the  calculation, 
the  feed  forward  controller  provided  the  Pulse  Width  Modulation 
(PWM)  command  through  the  PWM  module  and  directly  controlled 
the  air  blower  speed  so  that  the  flow  rate  at  the  cathode  inlet  could 
be  controlled. 
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Fig.  3.  Feed  forward  control  algorithm  for  PEM  fuel  cell  air  supply  control. 


Fig.  4.  Feed  forward  control  algorithm  for  PEM  fuel  cell  hydrogen  supply  control. 


4.2.  Hydrogen  flow  control 

Using  the  same  principle  as  the  air  supply  control,  the  hydrogen 
flow  control  was  achieved  using  a  simple  feed  forward  controller, 
as  shown  in  Fig.  4. 

The  feed  forward  coefficient  can  be  calculated  as  a  function  of 
the  load  following  using  Eq.  (33). 

Using  the  feed  forward  controller  for  cathode  and  anode 
flow  control,  only  one  sensor  is  required  to  complete  both 
tasks. 


4.3.  Simulation  of  PEM  fuel  cell  stack  system 

After  the  system  input  was  applied  (air  and  hydrogen  flows),  the 
simulation  results  of  the  PEM  fuel  cell  stack  system  with  integrated 
controller  were  available  for  analysis. 

The  parameters  used  in  this  simulation  were  based  on  a  large 
PEM  fuel  cell  stack  (40  kW  fuel  cell  stack)  to  more  clearly  interpret 
the  results. 

Fig.  5  shows  the  PEM  fuel  cell  polarization  characteristic  curve 
from  the  simulation  results.  This  curve  captures  the  real  PEM  fuel 
cell  polarization  behavior.  The  voltage  curve  decreases  as  the  fuel 
cell  current  density  increases  due  to  the  three  major  voltage  over¬ 
potentials.  The  power  curve  is  the  product  of  the  fuel  cell  current 
and  the  voltage.  Therefore,  system  power  first  increases  as  the  cur¬ 
rent  density  increases;  however,  when  the  current  density  reaches 
a  specific  value,  it  begins  to  decrease.  Eventually  the  system  power 
reaches  zero  when  the  output  voltage  reaches  zero. 

A  serial  step  input  signal  (Fig.  6)  was  employed  as  the  system 
input  to  test  the  fuel  cell  stack  system  model.  With  this  input,  the 
fuel  cell  voltage  response  to  the  current  demand  and  the  fuel  cell 
transient  could  be  investigated. 


Fig.  6.  Serial  step  load  input. 


The  fuel  cell  output  voltage  response  to  the  step  input  is 
shown  in  Fig.  7.  As  the  current  input  increases,  the  output  voltage 
decreases. 

The  PEM  fuel  cell  voltage  transient  response  is  shown  in  Fig.  8. 
This  figure  shows  the  fuel  cell  voltage  response  during  the  external 
load  change.  It  can  be  seen  from  this  figure  that  when  the  input 
data  changed  from  100  A  to  130  A,  the  voltage  did  not  reduce  to  a 
steady  level.  Instead,  it  dropped  to  a  lower  value  and  recovered  to 
a  steady  level  after  a  period  of  time.  This  phenomenon  was  caused 
by  the  slow  fuel  cell  transient  response. 

Fig.  9  shows  the  fuel  cell  output  power  response  to  the  cur¬ 
rent  step  input  increased.  However,  the  increasing  rate  tended  to 
be  slower  than  before.  This  is  because  as  the  demanded  current 
increases,  the  output  voltage  continues  to  decrease  due  to  the  fuel 
cell  polarization  behavior. 

Fig.  10  shows  the  relative  humidity  change  due  to  the  step  cur¬ 
rent  input.  Relative  humidity  is  an  important  parameter,  which 
directly  affects  the  fuel  cell  output  voltage.  In  the  fuel  cell  anode, 


Fig.  5.  PEM  fuel  cell  polarization  characteristic. 


Fig.  7.  Voltage  response  to  step  input. 
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the  RH  decreases  as  the  current  density  increases.  This  is  because 
the  electro-osmotic  drag  force  becomes  higher  than  the  diffusion 
drag  as  the  current  density  increases.  The  water  tends  to  move  from 
the  anode  to  the  cathode.  Also,  water  is  produced  during  the  fuel 
cell  operation;  therefore,  the  RH  in  the  cathode  is  always  equal  to 
one.  This  curve  can  be  used  for  humidification  control  of  the  fuel 
cell  operation. 

4.4.  Power  management  ofPEMfuel  cell-battery  power  system 
[28] 

In  this  study,  the  power  management  algorithm  is  proposed 
based  on  the  following  rules. 

When  the  load-required  power  is  lower  than  the  fuel 
cell-battery  system  rated  power,  which  is  100  W  in  this  study,  the 
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Fig.  10.  Relative  humidity  change  due  to  step  current  input. 
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Fig.  11.  Load-required  power. 


fuel  cell  will  be  employed  to  supply  sufficient  power  to  the  load.  The 
battery  will  only  be  used  to  supply  power  to  the  auxiliary  systems. 
If  the  battery  state  of  charge  (SOC)  is  lower  than  its  set  lower  bound, 
the  fuel  cell  will  start  to  charge  the  battery  until  the  SOC  reaches  its 
set  upper  bound.  During  the  charging  time,  the  auxiliary  systems 
will  be  powered  directly  by  the  fuel  cell. 

When  the  load-required  power  is  higher  than  the  system  rated 
power  but  lower  than  150%  of  the  system  rated  power,  the  battery 
will  be  controlled  to  supply  power  simultaneously  to  the  load  and 
the  fuel  cell.  The  battery  will  not  be  charged  whether  the  SOC  is 
lower  than  the  set  lower  bound  or  not.  Auxiliary  systems  will  be 
powered  by  the  battery  at  all  times. 

When  the  load-required  power  is  higher  than  150%  of  the  fuel 
cell-battery  power  system  rated  power,  the  system  will  automat¬ 
ically  shut  down  the  connection  to  the  load  to  prevent  damage  to 
the  load  device  and  the  fuel  cell-battery  power  system. 

4.4.1.  Simulation  conditions  of  power  management  system 

All  the  mathematical  models  introduced  above  were  imple¬ 
mented  in  the  Simulink  environment. 

To  capture  the  battery  charging  behavior  in  a  short  simulation 
time,  a  1 2  V,  0.1  Ah  lithium-ion  battery  was  used.  According  to  Kato 
in  [29],  the  ideal  SOC  of  a  lithium-ion  battery  should  be  maintained 
at  approximately  60%  to  extend  its  life.  Therefore,  in  this  study  we 
chose  the  SOC  upper  bound  to  be  60%,  the  lower  bound  to  be  58% 
and  the  initial  SOC  was  specified  at  60%. 

The  auxiliary-consumed  power  was  the  sum  of  the  power  con¬ 
sumed  by  the  blower,  the  solenoid  valves  and  the  microcontroller. 
The  power  consumed  by  the  microcontroller  was  assumed  to  be 
5  W  as  long  as  the  operation  continued. 

4.4.2.  Power  management  system  simulation  results  and 
discussion 

In  this  simulation,  a  series  of  step  signals  were  used  to  repre¬ 
sent  the  load-required  power,  which  is  the  input  to  the  system, 
as  shown  in  Fig.  11.  These  step  signals  can  be  used  to  represent 
most  cases  in  practical  applications.  At  time  0  s,  the  load  began  to 
require  power  from  the  power  system  at  10  W.  At  5  s,  the  required 
power  increased  to  60  W.  At  10  s,  the  required  power  reached  the 
power  system  rated  power  of  100W.  The  first  overload  power, 
1 40  W,  began  at  1 5  s ;  however,  this  was  within  the  system’s  afford¬ 
able  power  range,  which  was  150%  of  the  rated  power.  At  20  s,  the 
required  power  became  stable  at  the  rated  power  of  100W  until 
30  s.  At  30  s,  the  second  overload  power  of  180  W  occurred,  which 
exceeded  the  150%  system  rated  power  limit. 

The  simulation  results  are  shown  in  Figs.  11-18. 

Fig.  12  shows  the  fuel  cell  required  current  calculated  by  the 
power  management  system  and  sent  to  the  fuel  cell  stack  sub¬ 
system  as  an  input.  The  fuel  cell  stack  sub-system  took  the  required 
current  as  the  input  and  generated  the  fuel  cell  voltage  and  power, 
as  shown  in  Fig.  13.  The  fuel  cell  voltage  was  regulated  by  a  DC/DC 
converter  to  a  constant  12  V,  and  the  fuel  cell  power  was  reduced  by 
1 5%  due  to  the  efficiency  of  the  DC/DC  converter.  The  regulated  DC 
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Required  Fuel  Cell  Current 


Fig.  12.  Required  fuel  cell  current. 


Fig.  16.  Battery  SOC. 


Fuel  Cell  Produced  Power  and  Voltage 


Regulated  Fuel  Cell  Power  and  Voltage 


Fig.  14.  Regulated  fuel  cell  power  and  voltage. 


power  and  voltage  is  shown  in  Fig.  14.  The  battery  power  is  shown 
in  Fig.  15.  The  positive  power  value  represents  the  power  supplied 
by  the  battery,  and  the  negative  values  indicate  that  the  battery 
was  in  the  charging  process.  According  to  the  battery  power  situa¬ 
tion,  the  SOC  is  shown  in  Fig.  16.  The  auxiliary  systems  consumed 
power  over  time  as  shown  in  Fig.  17.  The  overall  output  power  vs. 
the  load-required  power  of  the  system  is  shown  in  Fig.  18. 

The  results  show  that  in  the  first  5  s,  the  load-required  power 
was  only  1 0  W  and  only  the  fuel  cell  provided  the  power  to  the  load. 
At  5  s,  the  load-required  power  increased  to  60  W.  The  required 
fuel  cell  power  increased;  thus,  the  battery  power  consumed  by 
the  blower  increased.  At  5.2  s  the  battery  SOC  reached  the  lower 
bound  of  58%,  and  the  battery  entered  the  charging  phase.  The  fuel 
cell  produced  power  increased  simultaneously  to  charge  the  bat¬ 


Battery  Power 


Auxiliary-Consumed  Power 


Fig.  17.  Auxiliary-consumed  power. 


tery  and  supply  the  auxiliary  systems.  At  8  s,  the  battery  SOC  was 
charged  to  60%;  therefore,  the  fuel  cell  ceased  charging  the  battery 
and  the  battery  began  to  supply  power  to  the  auxiliary  systems 
again.  At  10  s,  the  load-required  power  increased  to  the  system 
rated  power  of  100W.  The  fuel  cell  could  still  sufficiently  supply 
the  load-required  power  and  charge  the  battery  during  this  pro¬ 
cess.  At  15  s,  the  load-required  power  increased  to  MOW,  which  is 
40%  higher  than  the  system-rated  power.  During  this  situation,  the 
fuel  cell  was  not  able  to  provide  sufficient  power  for  the  load;  there¬ 
fore  the  power  management  controls  required  the  battery  to  supply 
power  to  the  load  simultaneously  with  the  fuel  cell.  At  1 8  s,  the  SOC 
was  below  the  lower  bound  of  58%,  and  the  battery  was  not  charged 
because  the  power  management  controls  prioritized  supplying  the 
load-required  load.  At  20  s,  the  load-required  power  was  reduced 
to  90  W,  which  was  lower  than  the  system-rated  power.  The  fuel 
cell  began  to  charge  the  battery  because  its  power  was  sufficient 
for  charging.  For  the  following  10  s,  the  load-required  power  was 
maintained  at  the  system  rated  power  to  ensure  that  the  system 
ran  in  a  stable  process.  At  30  s,  the  load  required  is  180W,  which 
is  more  than  150%  of  the  power  system-rated  power.  Therefore,  to 
protect  both  the  load  device  and  the  system,  the  overload  emer¬ 
gency  shut  down  process  was  launched.  The  power  system  ceased 
providing  power  to  the  load.  The  fuel  cell  only  provided  power  to 
the  battery  if  it  required  charging,  and  the  only  power-consuming 
device  in  use  was  the  microcontroller. 


Load-Required  Power  vs.  System  Output  Power 
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Fig.  15.  Battery  power. 


Fig.  18.  Load-required  power  vs.  system  produced  power. 
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5.  Portable  PEM  fuel  cell-battery  system  with  an  embedded 
control  system 

5.1.  System  explanation 

In  this  study,  a  100W  portable  PEM  fuel  cell-battery  system 
was  built  to  validate  the  designed  controller  performance  and  to 
investigate  the  PEM  fuel  cell  performance. 

A  portable  PEM  fuel  cell-battery  power  system  diagram  is 
shown  in  Fig.  19.  In  this  system,  a  microcontroller  was  employed 
as  the  ‘brain’  of  the  system.  The  microcontroller  used  the  demand 
current  as  an  input  and  used  a  power  management  algorithm  to 
calculate  the  power  distribution.  The  PEM  fuel  cell  stack  was  con¬ 
trolled  by  the  blower  and  the  hydrogen  control  valve.  The  battery 
charge  or  discharge  was  controlled  using  a  switch.  Temperature 
sensors  were  placed  within  the  system,  and  a  microcontroller 
monitored  the  system  temperature  through  these  sensors  and 
sent  a  control  signal  to  a  cooling  fan  to  control  the  system 
temperature. 


5.2.  Validation  of  the  PEM  fuel  cell  control  system  simulation 
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Fig.  20.  Step  current  input  for  simulation  and  experiment. 


such  as  rapid  temperature  and  relative  humidity  changes,  we  chose 
a  low  current  density  zone  to  test  the  fuel  cell  stack. 

Fig.  21  shows  the  voltage  output  of  the  system  for  the  sim¬ 
ulation  and  the  experiment.  The  simulation  results  follow  the 


In  this  study,  the  100W  PEM  fuel  cell-battery  power  system 
with  an  embedded  control  system  was  used  to  validate  the  simu¬ 
lation  results.  The  fuel  cell  stack  used  here  was  a  100  W  PEM  fuel 
cell  stack  with  10  single  cells  and  a  cross-sectional  area  of  50  cm2. 

The  simulation  environment  was  set  up  as  the  experimental 
environment  with  the  following  parameters: 

i.  room  temperature, 

ii.  an  initially  well  humidified  stack  was  established  by  injecting 
liquid  water,  and 

iii.  the  same  inlet  initial  pressure  and  ambient  outlet  pressure. 

The  system  inputs  for  the  simulation  and  the  experiment  are 
shown  in  Fig.  20.  A  step  current  signal  was  employed.  To  obtain 
more  accurate  and  clear  results  by  minimizing  the  uncertainties, 


Fig.  21.  Voltage  output  for  simulation  and  experiment. 
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Fig.  22.  Power  output  for  simulation  and  experiment. 


Fig.  23.  Step  current  demand  input. 


trend  of  the  experimental  results.  They  correspond  well  for  the 
first  current  (1A)  because  in  the  initial  state,  all  conditions  in 
the  simulation  and  the  experiment  were  the  same  as  defined 
at  initialization.  As  the  demand  current  increased,  the  experi¬ 
mental  results  became  lower  than  the  simulation  results.  The 
gap  between  the  simulation  results  and  the  experimental  result 
increased.  This  phenomenon  was  caused  by  two  issues.  (1)  As  the 
experiment  continued,  the  relative  humidity  inside  fuel  cell  stack 
decreased;  thus,  the  output  voltage  decreased  more  than  the  sim¬ 
ulation  results.  The  decreasing  relative  humidity  also  caused  the 
fuel  cell  operation  response  time  to  slow.  (2)  In  the  experiment, 
a  12  V  DC  2  A  blower  was  used  as  the  air  supply  tool.  The  max¬ 
imum  output  pressure  of  this  blower  was  0.5  psi,  which  meets 
the  minimum  requirement  for  this  fuel  cell  stack.  However,  the 
fuel  cell  stack  had  a  long,  thin  channel,  which  caused  a  signif¬ 
icant  pressure  drop.  Therefore,  the  actual  cathode  pressure  was 
significantly  lower  than  the  blower  output  pressure.  When  the 
demand  current  was  1  A,  the  input  flow  rate  of  the  fuel  cell  stack 
was  sufficient;  however,  when  the  demand  current  increased,  the 
fuel  cell  did  not  have  sufficient  oxygen  as  a  reactant.  Thus,  the 
experimental  output  voltage  decreased  faster  than  the  simula¬ 
tion  results.  This  problem  can  be  solved  using  a  high-pressure 
blower;  however,  the  larger  blower  will  consume  more  power  than 
the  current  blower.  Because  the  blower  power  was  provided  by 
the  fuel  cell  stack  output  power,  it  had  a  theoretical  power  limit 
of  100W.  Currently,  a  12  V  blower  that  can  provide  up  to  5psig 
of  pressure  and  consume  less  than  100W  of  power  cannot  be 
found. 

Fig.  22  shows  the  power  output  for  the  simulation  and  the 
experiment.  It  follows  the  trend  of  the  system  voltage  output.  The 
increasing  speed  of  the  experimental  results  was  becoming  slower 
than  the  simulation  results.  The  difference  was  also  caused  by  air 
starvation  and  low  relative  humidity. 

5.3.  Portable  PEMfuel  cell-battery  power  system  performance 

The  purpose  of  building  a  100  W  portable  PEM  fuel  cell-battery 
power  system  was  to  use  the  PEM  fuel  cell  stack  as  the  main  power 
source  to  provide  continuous  regulated  DC  power  to  electronics. 
Therefore,  in  the  testing  of  the  power  system,  we  tested  the  reg¬ 
ulated  DC  power  output  from  the  PEM  fuel  cell-battery  power 
system.  The  expected  result  was  a  constant  12  VDC  output  voltage, 
which  is  not  affected  by  the  current  demand  like  the  fuel  cell  stack 
system.  Therefore,  this  result  did  not  require  simulation  results  for 
validation. 

A  step  current  signal  (shown  in  Fig.  23)  was  used  in  this  exper¬ 
imental  study  as  the  system  input. 

Fig.  24  shows  the  system  voltage  output  of  the  power  system. 
For  the  first  1 80  s,  the  system  output  voltage  was  between  12V  and 


11.8  V,  which  is  acceptable  for  most  electronics.  However,  begin¬ 
ning  at  180  s,  the  system  output  became  unstable.  The  instability 
was  caused  by  the  properties  of  the  DC/DC  converter.  The  DC/DC 
converter  used  in  this  work  accepts  an  input  voltage  from  6  to 
19  V  DC,  and  the  efficiency  is  approximately  80%,  when  the  input 
voltage  is  between  9  and  15  V.  When  the  current  demand  at  the 
DC/DC  converter  output  was  2  A,  the  actual  current  demand  of  the 
fuel  cell  stack  was  5  A  or  more  because  the  fuel  cell  stack  output 
voltage  at  this  time  is  approximately  7  V.  At  this  point,  the  DC/DC 
converter  efficiency  becomes  non-linear,  which  can  lead  to  an  effi¬ 
ciency  less  than  30%.  The  higher  the  current  demand  on  the  fuel 
cell  stack,  the  lower  the  fuel  cell  stack  voltage  becomes.  When 
the  fuel  cell  stack  output  voltage  was  less  than  6  V,  the  DC/DC 
converter  ceased  running.  No  commercial  DC/DC  converter  that 
could  accept  an  input  voltage  less  than  6  V  was  found,  as  physi¬ 
cally  converting  such  a  low  voltage  to  a  high  voltage  will  cause  a 
huge  power  loss.  Therefore,  the  solution  is  either  to  build  a  DC/DC 
converter  or  use  a  fuel  cell  stack  that  can  provide  a  higher  output 
voltage. 

Fig.  25  shows  the  fuel  cell-battery  power  system  output  power. 
Due  to  the  limitation  of  the  low  voltage  output  associated  with  the 
fuel  cell  stack,  the  system  output  cannot  be  higher  than  1 8  W.  How¬ 
ever,  as  a  demonstration  project,  these  results  show  that  the  system 
can  provide  a  continuous  12  V  DC  power  supply  to  electronics  as  a 
power  source. 

From  the  results  in  Section  6,  it  is  clear  that  there  are  some 
physical  limitations  of  the  hardware  that  prevent  the  PEM  fuel 
cell-battery  power  system  from  achieving  an  improved  per¬ 
formance.  A  high-pressure  blower  or  air  pump  is  required  to 
prevent  the  fuel  cell  stack  from  air  starvation  when  the  demand 
current  increases.  To  power  the  high-pressure  blower  or  air 
pump  and  to  increase  the  output  voltage  to  fit  the  work- 
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ing  voltage  range  of  the  DC/DC  converter,  a  larger  fuel  cell 
stack  with  a  higher  output  voltage  is  recommended  for  further 
study. 

Also,  it  is  clear  from  the  experimental  results  that  the  fuel  cell 
output  voltage  response  tends  to  be  slower  as  the  relative  humidity 
decreases.  A  physical  humidification  system  is  required  to  improve 
the  fuel  cell  system  performance. 

6.  Conclusions 

In  this  study,  a  general  methodology  of  modeling,  control  and 
building  a  proton  exchange  membrane  fuel  cell-battery  power  sys¬ 
tem  was  developed.  A  set  of  PEM  fuel  cell-battery  power  system 
models  was  introduced.  This  model  can  be  used  to  address  the 
transience  of  the  PEM  fuel  cell-battery  power  system  and  can  be 
used  as  a  platform  for  system  controller  development.  A  parametric 
study  of  the  influences  of  the  temperature,  the  relative  humidity, 
the  inlet  pressure,  and  the  oxygen  excess  ratio  on  the  fuel  cell  per¬ 
formance  was  conducted  based  on  this  model.  A  set  of  PEM  fuel 
cell-battery  power  system  controllers  was  developed.  The  results 
show  that  the  air  supply  control  system  could  deliver  sufficient 
reactants  for  the  fuel  cell  stack;  the  humidification  control  system 
was  able  to  increase  the  fuel  cell  system  performance  by  main¬ 
taining  the  fuel  cell  stack  relative  humidity;  and  the  system  power 
management  system  was  able  to  accurately  distribute  the  demand 
power  between  the  fuel  cell  and  the  battery  and  maintain  the  sys¬ 
tem  components  life  and  emergency  shutdown.  A  physical  PEM 
fuel  cell-power  system  with  embedded  control  system  was  built 
to  validate  the  simulation  results  and  to  demonstrate  the  system. 


The  validation  results  show  that  the  simulation  results  correspond 
to  the  experiment  results. 
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